Zeeman Relaxation of Cold Atomic Iron and Nickel in Collisions with ^He 
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We have measured the ratio 7 of the diffusion cross-section to the angular momentum reorientation 
cross-section in the colliding Fe-'^He and Ni-^He systems. Nickel (Ni) and iron (Fe) atoms are 
introduced via laser ablation into a cryogenically cooled experimental cell containing cold (< 1 K) 
''He buffer gas. Elastic collisions rapidly cool the translational temperature of the ablated atoms to 
the ^He temperature. 7 is extracted by measuring the decays of the atomic Zeeman sublevels. For 
our experimental conditions, thermal energy is comparable to the Zeeman splitting. As a result, 
thermal excitations between Zeeman sublevels significantly impact the observed decay. To determine 
7 accurately, we introduce a model of Zeeman state dynamics that includes thermal excitations. We 
find 7Ni-3Hc = 5x 10^ and 7fo_3hc < 3x 10^ at 0.75 K in a 0.8 T magnetic field. These measurements 
are interpreted in the context of submerged shell suppression of spin relaxation as studied previously 
in transition metals and rare earth atoms. [H-Q- 
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I. INTRODUCTION 

Cooling and trapping atoms at cold and ultracold tem- 
peratures has led to unprecedented control of the exter- 
nal and internal atomic degrees of freedom. Such con- 
trol has aided in the advance of quantum information, 
precision measurement and atomic clocks, and has set 
the stage for quantum simulation of condensed matter 
systems. Successful evaporative cooling of trapped, di- 
lute atomic gases requires a thorough understanding of 
elastic collisions that thermalize the sample and inelas- 
tic processes that cause heating and atom loss. Efficient 
evaporation is possible only for species with a high ra- 
tio of elastic to inelastic collisions 7. For atoms with 
isotropic interactions, such as the alkalis, 7 is typically 
large. Alkalis can also be conveniently laser cooled due 
to their simple electronic structure. For these reasons, 
alkalis have been the atoms of choice for the majority of 
ultracold atomic physics experiments. 

Because of the proliferation of applications for ultra- 
cold atoms, there is a natural desire to develop techniques 
to expand the class of atomic species that can be cooled 
and trapped to those with rich electronic structures. Re- 
cently, alkaline earth metal atoms and atoms with similar 
electron structure have been cooled to quantum degener- 
acy using a combination of laser cooling and evaporation 
out of an optical dipole trap. Examples include *°Ca [1], 
^■^Sr [5|,|6|], and multiple ytterbium isotopes 0-[ll|. There 
is also great interest in systems with large anisotropic 
magnetic dipole-dipole interactions, which have been pre- 
dicted to have novel quantum phases [l^ . have poten- 
tial use in quantum computing |13l|. a nd demonstrate 
geometry-dependent BEC stability [Tsj . Although 
metastable states of alkali earth atoms are one candi- 
date to demonstrate the above effects [l^ , the measured 
inelastic loss rates in trapped samples are high [iMl, 



making evaporative cooling difficult. Identifying viable 
alternatives with even larger magnetic moments requires 
the development of more general methods of trapping 
and cooling exotic species. 

One technique capable of trapping exotic species is 
buffer-gas cooling ^20.]. In buffer-gas cooling, elastic colli- 
sions with a cold buffer gas, typically helium, are used to 
cool the atom or molecule of interest to low enough tem- 
peratures that the species may be magnetically trapped. 
Traditional evaporative cooling techniques may then be 
used to cool the species to the regime of ultracold tem- 
peratures. The flexibility of buffer-gas cooling has been 
demonstrated by the successful trapping of many non- 
alkali atomic and molecular species [3, [2l| - |27| . Further- 
more, a buffer-gas cooled Bose-Einstein condensate of 
metastable helium {'^He*) was recently reported [28| , 
demonstrating the viability of the technique to create 
ultracold gases. Buffer-gas loading of magnetic traps re- 
quires about 50 collisions in order for the target species 
to fully thermalize with the buffer gas. However, it takes 
additional time for the atoms to diffuse toward the center 
of the trap and for the buffer gas to be removed. Thus 
buffer-gas loading of magnetic traps has required that 
the orientation of the species' magnetic moment must be 
preserved for at least 10"' collisions 

To first order, angular momentum reorientation in such 
collisions is dominated by the interaction between the or- 
bital angular momentum L of the species and the angular 
momentum i of the collision [s^l- One therefore expects 
weak reorientation for species having L = 0, and rela- 
tively strong reorientation for species having L ^ 0. For 
instance, oxygen (L = 1) reorients its angular momentum 
in almost every collision with low-temperature He [30| , 
while atoms such as potassium or copper (L — 0) have 
been shown to survive more than 10^ collisions without 
angular momentum reorientation (22l . [sij . Recent work 
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[32l |33| . however, has demonstrated the existence of a 
class of atoms, dubbed the "submerged shell" atoms, m 
which the L ^ Q valence shell is protected from angu- 
lar reorientation by outer-lying filled shells of electrons. 
In these species, the anisotropy of the interaction with 
helium is highly suppressed and 7 is orders of magni- 
tude larger than typically found in non-S-state atoms. 
Experimental work with rare earth atoms Q and the 
group 3-5 transition metals [l[ has shown the ratio 7 
of diffusion cross-section to angular momentum reorien- 
tation cross-section to be between 10^ and 10^ in these 
"submerged shell" atoms. Similar suppression of Zee- 
man relaxation has been observed in collisions between 
helium and species due to their spherically sym- 

metric electron-density distribution [34]. Suppression of 
fine-structure-changing collisions has also been observed 
in collisions between "submerged shell" transition metal 
titanium with helium (ssj . 

Previously only a few transition metals were studied [l| 
and each species had a small magnetic moment (< 1.32 
bohr magnetons). Because species with large magnetic 
moments have stronger dipole-dipole interactions and 
are easier to trap, characterizing transition metals with 
larger magnetic moments is desirable. Our apparatus is 
equipped with a large cryogenic valve that can be rapidly 
opened to remove the buffer gas after initial cooling, al- 
lowing us to thermally isolate and evaporatively cool a 
trapped sample. Thus we are well positioned to study 
species with 7 in the range of previously studied "sub- 
merged shell" atoms. 

Our goal was to investigate the suppression of reorien- 
tation for transition metals with large magnetic moments 
in their collisions with low-temperature '^He. Specifically, 
we study 7 for the Ni-'^He and Fe-^He systems because 
nickel and iron have strong spectroscopic lines accessi- 
ble to our doubled dye laser system^ . We find that these 
atoms have 7 < 5 x 10^. As with previously studied tran- 
sition metals, their reorientation is more rapid than rare 
earth "submerged shell" species and is, unfortunately, 
too rapid for further study in our apparatus. However, 
we find that nickel (7 ~ 5 x 10'^) still demonstrates sig- 
nificant reorientation suppression compared to strongly 
anisotropic species. 

When the thermal energy of the atomic sample is com- 
parable to the energy separation between adjacent Zee- 
man states, thermal excitations have a non- negligible ef- 
fect on Zeeman state dynamics. This manifests itself as 
a decrease of the observed decay rate compared to the 
angular momentum reorientation decay rate. As a re- 
sult, a naive model that neglects these effects will lead 
to an overestimation of 7. We introduce a method for 
extracting 7 from observed Zeeman state decay by in- 



^ Cobalt was not included in this study as the large hyperfine 
structure of the atom precluded spectroscopic identification of 
the atom's individual Zeeman states. 



Atom 


Configuration 


Term 


Moment {^b) 


Fe 


[Ar].3d^4s2 




6.005 


Ni 


[Ar].3d*.4s2 




5.002 



TABLE I: Electronic configurations of iron and nickel. 

eluding thermal excitations in our Zeeman state dynam- 
ics model. We also discuss the impact of our method 
on the interpretation of previous measurements of 7 for 
transition metals and rare earth atoms. 



II. EXPERIMENTAL METHODS 

We are interested in the ratio between the atom-He dif- 
fusion cross-section ad and the angular momentum reori- 
entation cross-section a^. We are specifically concerned 
with reorientation from the most low-field seeking trap 
state, mj = J, to any other Zeeman state. In the pres- 
ence of a confining magnetic field, these other states will 
leave the magnetic trap on a time scale exponentially 
faster than the lifetime of the mj = J state. We there- 
fore write 

<^R= X! '^J^m'j, (1) 

7-^. (2) 

The diffusion cross-section may be measured by observ- 
ing the diffusion of the atoms to the wall when no mag- 
netic field is present. To measure cr_R, we apply an ap- 
proximately uniform magnetic field to separate the var- 
ious mj sublevels of the atom, then measure the time 
constant for loss of the mj = J population as a function 
of the zero- field diffusion time. Extraction of the cross- 
section ratio from measurements of the atom decay time 
constants is discussed in ijllll 

A. Cryogenic Apparatus 

A schematic representation of our apparatus (not 
drawn to scale) is shown in Fig. [TJ The body of the 
experimental cell is machined from a GIO tube 7.6 cm in 
diameter and 30 cm in length. It is thermally anchored 
to the mixing chamber of a dilution refrigerator by four 
half-inch diameter oxygen free high purity copper braids. 
Thermal conductivity along the length of the cell is pro- 
vided by ^ 1,000 0.25 mm diameter copper wires running 
vertically along the outer GIO cell wall. Each wire is elec- 
trically insulated from the rest to prevent eddy current 
heating as we ramp magnetic fields. Base temperature 
of the refrigerator is 30 mK resulting in a cell top tem- 
perature of ~ 100 mK and cell bottom temperature of ~ 
170 mK. A 6.6 cm sapphire window forms the bottom of 
the cell. Sapphire passes the UV frequencies needed for 
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FIG. 1: (color online) Schematic drawing of the cryogenic 
apparatus. The magnet can operate in Helmholtz (shown) or 
anti-Helmholtz configurations. 



spectroscopy and has higher thermal conductivity than 
fused silica. The window is epoxied into the GIO body, 
forming a vacuum seal with the cell wall. 

The cell body is housed in a cylindrical vacuum cham- 
ber with 1 mm radial clearance. A superconducting mag- 
net fits tightly around the vacuum chamber [36j. Clear- 
ances are made as tight as possible to enable maximum 
fields at the cell wall. The room temperature connec- 
tions to the magnet coils can be wired to produce a 
Helmholtz field, providing a uniform field within the cell, 
or an anti-Helmholtz field, providing a 4 T deep spherical 
quadrupole trap within the cell. This work required only 
Helmholtz fields. 

Buffer gas is introduced into the cell through a fill line. 
The fill line is connected to a vacuum chamber containing 
a charcoal sorb filled with ■^He. The vacuum chamber is 
thermally anchored to a 4K helium bath. For each buffer- 
gas load, we heat the sorb with a resistive heater to 10 
K. This drives helium off the sorb, through the fill line, 
and into the cell. Adjusting the duration and power of 
the heater pulse allows us to vary the amount of loaded 
buffer gas in a predictable manner. Before each decay 
lifetime measurement is performed, the cell is heated to 
~ 350 mK. This drives the buffer gas off the walls in 
preparation for introducing the atomic species. 



B. Spectroscopy 

Fe and Ni atoms are produced via ablation of metallic 
targets mounted inside the cell. Ablation is performed 
with a 10 ns pulse from a doubled YAG laser operating 
at 532 nm. Both atom density and temperature increase 
with ablation power. We used pulse energies ~ 15 mJ 
to reach densities that yielded adequate signals. Un- 
fortunately, these powers also resulted in temperatures 
at which thermal excitations between mj states signifi- 
cantly contributed to the observed decay rates. 

We probe atomic density, lifetime, and temperature 
via balanced absorption spectroscopy on the ^Z?4 ^F^, 
transition at 248 nm in Fe and the — ^G^ transition 
at 232 nm in Ni. The light in both cases was produced 
from a dye laser doubled in a resonant cavity containing 
a BBO crystal. 

Optical access into the cryogenic apparatus is limited 
to a single port through the bottom. Beam steering op- 
tics mounted to the bottom of the dewar direct the laser 
into the cell. The beam retroreflects from a mirror at 
the top of the cell and the exiting light is detected on a 
photomultipler tube (Hamamatsu H6780-04 [st}). 



C. Measurement of the diffusion lifetime 

The lifetime of the atoms at zero-field is set by their 
diffusion through the buffer gas to the cell wall. The dif- 
fusion lifetime is proportional to the buffer-gas density 
nb in the cell (see equation U]) . A measurement of is 
therefore a relative measurement of ni, . 

To make the lifetime measurement the frequency of the 
laser is tuned to the atomic resonance of interest. The 
laser frequency is scanned repeatedly over the absorption 
feature as atoms are introduced into the cell. The number 
of absorbing atoms in the probe beam is measured by 
integrating the spectrum over the atomic line. Drifts 
in the other sources of loss in the optical path can be 
mitigated by subtracting the signal at a dark portion of 
the spectrum from this integral. We fit the decay in this 
integrated optical depth to a single exponential decay 
function to determine the diffusion lifetime r^. Because 
the laser scan has a maximum bandwidth of ^ 30 Hz, 
decays faster than this are measured by parking the laser 
frequency at the absorption peak. However, this latter 
method is subject to low frequency noise and drift. 



D. Measurement of the mj = J lifetime 

To measure the lifetime tj of the mj ~ J state, we 
must resolve an absorption peak from the mj — J ground 
state. This is accomplished by turning on a homogeneous 
magnetic field, thereby splitting the atom's Zeeman sub- 
levels. The lines are also broadened, due to field inhomo- 
geneity within the magnet. 
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Ni Amj = 1 transition, 0.8 T field 



58Ni nij = 3 
''"Ni mj = 4 



s^Ni mj = -4 
eoNi mj = -3 
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FIG. 2: Nickel optical depth vs. frequency (arbitrary zero) in 
a homogeneous (Helmholtz) field. Each isotope has 9 lines 
corresponding to the 9 mj states. The Zeeman and isotope 
shifts are roughly equal at 0.8 T fields, causing lines of differ- 
ent isotopes to overlap. Measurements of the mj — J state 
lifetime are performed by parking on the mj = 4 transition 
peak and measuring the optical depth vs. time. 



Consider a transition from a mj ground state with 
Lande factor gj to a m j + Am j excited state with Lande 
factor g'j. The field induced frequency shift Avz is 



^'^z = (gjAmj + {gj - gj)mj 



flBB 



(3) 



where h is Planck's constant. Selection rules require 
Amj = 0, ±1. The first term in is generally 
much larger than the second, yielding three manifolds 
of equally spaced peaks, each corresponding to a ground 
state Zeeman sublevel. 

The mj — J state is identified by tuning to the 
Amj = 1 set of peaks as shown in Fig. [51^ Nickel's spec- 
trum in a magnetic field is relatively simple because the 
most common isotopes, ^^Nij^^Ni, and ^^Ni, have no nu- 
clear spin and hyperfine effects are absent. Each isotope 
splits into 9 lines, corresponding to the mj sublevels. 
The Zeeman splitting at ~ 0.8 T is approximately 
equal to the isotope shift between ^^Ni and ^"Ni. As a 
result, transitions from the ^^Ni mj state overlap with 
^"Ni transitions from the mj + 1 state. Only the ^^Ni 
mj — J and ^°Ni mj = — J states do not experience any 
overlap. We measure tj by tuning the laser frequency 
to the fully resolved ^^Ni mj = 4 absorption peak and 
observing the decay of the absorption signal. The most 
abundant iron isotopes (^^Fe and ^''Fe) are also J = 4 
species without hyperfine structure, so the spectroscopic 
methods outlined above also apply to Fe. We choose to 
operate at 0.8 T with Ni and 1.0 T with Fe, as these 



^ Clebsch-Gordan coefficients suppress the mj 
Amj = —1 and Amj = manifolds. 



J peak in the 



are the highest fields at which the spectra are easily in- 
terpreted and at which the absorption lines are not too 
broadened to achieve a good signal-to-noise ratio. 

Measuring tj was particularly difficult for iron because 
the observed lifetimes were very short at all buffer-gas 
densities. Immediately after ablation we typically see an 
absorption signal that decays in 2-3 milliseconds. This 
signal occurs at all buffer-gas densities and is present 
when the laser is parked near, but not necessarily on, a 
Zeeman level resonance peak. We associate this transient 
signal with the decay of higher order diffusion modes. 
Because it is unrelated to the momentum reorientation 
we are trying to measure, we typically ignore data taken 
in the first 10 ms after ablation. For example, Ni rj were 
found using data taken between 10 ms - 50 ms after 
ablation. However, due to Fe's rapid Zeeman relaxation, 
the absorption signal was too small to be measured after 
10 ms. We therefore measured rj for Fe starting only 5 
ms after ablation. 

Using the above procedure we can also study the dy- 
namics of the mj — —J most high-field seeking state. 
After several Zeeman relaxation times have elapsed, ther- 
mal equilibrium is established between the Zeeman and 
kinetic degrees of freedom and the mj = — J state decays 
via diffusion. We measure r_j at each buffer-gas density 
by observing mj = — J decay after thermal equilibrium 
has been established. We then compare t_j with t^. Un- 
der these conditions we expect t_j to equal r^. Fig. [3] 
shows an example of mj = J and mj = — J state decay 
with fits to a single exponential lifetime. 



E. Temperature Measurement 

At atom temperatures greater than or comparable to 
the atomic Zeeman splitting, 670 mK at 0.8 T for Ni, 
thermal excitations cause the observed loss rate of the 
mj = J state to differ from the Zeeman relaxation rate. 
In order to know if these excitations can be ignored, we 
determine atom temperature by measuring the broaden- 
ing of the zero-field spectrum. Fig. 2] shows the zero-field 
spectrum of the a^F4^ V^G^ transition of Ni at 232 nm, 
taken 50 ms after ablation. The optical detuning is cal- 
ibrated using a Fabry-Perot cavity. The observed atom 
density is 3 x 10* cm^"^, corresponding to 3 x lO^'' 
Ni atoms in the cell. The temperature of the atoms is 
determined by fitting to a Voigt profile [38j . 

Fig.[5]shows the temperature of the Ni atoms as a func- 
tion of time for three buffer-gas densities. Thermal exci- 
tations have the greatest impact at high buffer-gas den- 
sities, where the Zeeman relaxation rate is much faster 
than the diffusion rate. In our analysis, we shall use the 
average measured Ni temperature at high densities, 750 
mK. 

The zero-field spectrum could potentially be broadened 
due to trapped fluxes in the superconducting magnet. 
These have previously been measured to be < 10 gauss 
[39| . yielding a systematic uncertainty of mK. 
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FIG. 3: (color online) Nickel m,j = J and mj — — J optical 
depth (integrated over the atomic line) vs. time. Upper plot 
shows Zeeman state evolution starting immediately after laser 
ablation. The mj — — J state experiences an initial increase 
in optical depth because mj ^ —J states are relaxing into 
the mj = —J state. Lifetime fits are performed over the 
indicated regions. Lower left(right) plot shows mj — J{ — J) 
lifetime fit. 
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FIG. 4: (color online) Nickel zero-field optical depth vs. fre- 
quency 50 ms after ablation. Frequency zero set to ^*Ni res- 
onance. Temperature and density are found by fitting to a 
Voigt profile. 
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FIG. 5: (color online) Nickel temperature vs. time. Tem- 
perature measurements were made at 3 different buffer-gas 
densities. The atoms cool only slightly over the time scale of 
our Tj measurements. The temperature is slightly higher at 
lower buffer-gas densities. 



III. LIFETIME MODELS 
A. Model in the near-zero temperature limit 

At high bufFer-gas densities and in the limit of near- 
zero temperature, atoms are lost from the mj — J state 
due to two primary mechanisms. First, atoms diffuse 
via elastic collisions until they reach the cell walls, where 
they stick. Second, atoms experience Zeeman relaxation. 
In this section we ignore collisions that repopulate the 
mj — J state once the atom has relaxed into a lower en- 
ergy state. We will include this finite temperature effect 
in Section HiTbI 

In a cylindrical cell of radius i?ccii and length Lcoii, the 
lifetime due to diffusion for atoms in the lowest order 
diffusion mode is 140 



Td 

32 



rib <yd 
vG ' 

.^2 



(4) 
(5) 



where joi — 2.40483 ... is the first zero of the Bessel 
function Jo(2:), is the buffer-gas density, and v — 
{SkBT/inrf/^ is the mean relative velocity of the col- 
liding atom-^He system with reduced mass fx. We can 
ensure that all atoms are in the lowest order diffusion 
mode by a waiting for a few multiples of Td before mea- 
suring the atom lifetime. 

Near zero temperature, the lifetime of atoms in the 
m,j = J state is the reciprocal sum of the diffusion life- 
time and the lifetime due to Zeeman relaxation: 



1 1 
Tj = I — + r 

Td VUBUb 



(6) 



Unfortunately we do not have an absolute calibration of 
rib. However, we can use (U) to eliminate rif, from ([5]). 



6 



Doing this, and replacing Ud/o'R with 7, we find 



Tj 



1 + r| 1-2 G/7 



(7) 



7 can be found by measuring tj as a function of (mea- 
sured at zero magnetic field), and fitting to ([7]). However, 
as we shall show in the next section, this model breaks 
down when /csT > qj^bB. 





"Amj 




a±i 


±1, 2 


1.0 


±3, 4 


0.2 


±5, 6 


0.04 


±7, 8 


0.008 



TABLE II: Relative rate coefficients for Zeeman rela^cation 
with a given Amj, based on [30I. |4^. 



B. Model at finite temperature 

The previous model assumes that once a mj = J atom 
experiences Zeeman relaxation it remains in a lower en- 
ergy Zeeman state forever; the possibility of excitation 
into the mj = J states was ignored. When the ther- 
mal energy ksT is much less than the magnetic interac- 
tion energy, this assumption is valid as collisions do not 
have sufficient energy to excite atoms into states with 
higher mj. However, for our experimental conditions 
ksT << gjfisB does not hold, and a non-negligible per- 
centage of collisions with the buffer gas have enough en- 
ergy to excite an atom into a higher energy Zeeman state. 
This effect slows the observed loss from the mj — J state, 
yielding an overestimate of 7 when the near-zero temper- 
ature model is used. Therefore, a correct determination 
of 7 from measured data must take thermal excitations 
into account. 

For the above reasons a Zeeman cascade model which 
includes diffusion, the dynamics of all Zeeman states, and 
the possibility of thermal excitation must be developed. 
In our model, the density n^j of each mj level with 
energy Em, at temperature T evolves according to 

rn'j<mj 



r_R ^ Ctm'jm., nmj GXp 
m'j>7n .7 



m'j'yrn j 



Tr ^ a™>,n™'^exp( — 

m'j<mj 



(8) 



where Yd = 1/Td and Tn = 1/tji = ajiribv are the diffu- 
sion and relaxation rates. 7 is related to these quantities 
hy J ^ Til TdGv'^. o-m' m , in each summation represents 
the coupling between mj and m'j Zeeman levels, subject 
to 



2J 



and 



(9) 



(10) 



The first term in is diffusion loss, the second is Zee- 
man relaxation into lower energy states, the third is ther- 
mal excitation into higher energy states, the fourth is 
Zeeman relaxation from higher energy states, and the 
fifth is thermal excitation from lower energy states. 

The dominant effect of thermal excitation on mj ~ J 
atoms is an alteration of their initial decay from the sim- 
ple exponential predicted by the near-zero temperature 
model. At early times {t ^ l/r^j, when all mj states 
are equally populated), thermal excitations cause the 
nij — J state to decay more slowly than the prediction 
of ([7]). The exact scaling will depend on T, B, and the 
exact form of a, 



m jTYij 



At later times (t > (2 J + 1)/Tb), 
the Zeeman states approach thermal equilibrium, and the 
lifetimes of all mj levels approach r^. 

Finding reliable values of am'jmj presents a challenge 
[4l|. Because the rates for iron and nickel are unknown, 
they must be estimated. The assumptions adopted sig- 
nificantly affect the predicted Zeeman state dynamics. 
Consequently, we have analyzed them using three sce- 
narios for transitions between Zeeman states: all tran- 
sitions between Zeeman states are equally allowed, only 
Amj = ±1 transitions are allowed, and an intermediate 
regime based on calculations for thulium and general 
atoms [S^, \^ . The relative values of a used in this third 
case are shown in Table [TTl 



IV. DETERMINATION OF 7 

We use the finite temperature model to fit for 7 using 
the following method: For each experimentally observed 
value of Td, we simulate the mj = J state decay using a 
guess value of 7, a temperature of 750 mK, a magnetic 
field of 0.8T, and the literature-based a values from Ta- 
ble ini We then fit the simulated decay of the mj — J 
state over the same time interval used to measure tj. Fi- 
nally, we perform a fit of the simulated values of tj 
to measured values to determine 7. 



A. 7 for Ni— ' He collisions 

The TOj — J state lifetimes are plotted with the 
nij = —J state lifetimes in Fig. [Bl As expected, the 
TO J = — J state lifetimes are approximately Td, whereas 
the TOj = J state lifetimes first increase, then decrease 
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FIG. 6: (color online) Nickel mj — J and mj = — J lifetimes 
in a 0.8 Tesla field vs. diffusion lifetime. The dashed hne was 
drawn with a slope = 1 to demonstrate that the mj — —J 
atoms leave the cell by diffusion as expected. The mj — J 
atoms decay quickly due Zeeman relaxation in the region of 
high Td- The best fit of the mj = J data to a finite tempera- 
ture Zeeman cascade simulation yields a value of 7 = 5 x 10"^ . 

with increasing r^. The clear differentiation between the 
TO J = J state and mj = — J state behaviors provides 
convincing evidence that we are measuring Zeeman re- 
laxation of the TOj — J state. 

Fig. |6] shows the best fit of the mj — J data to the 
finite-temperature model, yielding 7 = 5 x lO''. This is 
nearly a factor of two smaller than the 7 = 9x10'^ ob- 
tained when fitting to the near-zero temperature model. 
Uncertainty in 7 results from three primary effects: de- 
viation of the data from the model, uncertainty in ther- 
mal excitation rates due to temperature uncertainty, and 
uncertainty in the assumed "selection rules" for Zeeman 
relaxation. 

Our data generally show deviations from the model 
larger than their statistical uncertainties. At low r^, we 
believe this is due to an unknown experimental artifact, 
as was also observed in our experiments with Cu-He and 
Ag-He [l^]. At high Td, deviation is caused by perform- 
ing measurements in a regime where loss is not purely 
exponential, as we begin to see the atomic states ap- 
proach their steady-state distributions. We account for 
the effect of this model deviation on our measurement of 
7 using an F test f43!|, yielding a parameter uncertainty 
with a 95 % confidence interval of x 10'^. 

Finally we consider the effects of the assumed selection 
rules for Zeeman relaxation. For the fit in Fig. [6] we used 
guidance from the literature in setting realistic selection 
rules. Since it is impossible to quantify the accuracy 
of this assumption, we systematically vary the selection 



Selection Rule 


Fit for 7 


From Literature 


5 X 


10^ 


All Equal 


7 X 




Am J — 1 only 


2 X 


10^ 



TABLE 111: Values of 7 extracted from fits of data to Zeeman 
cascade simulations under various selection rules assumptions. 
For these fits, we used T = 0.75 K. 
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10 20 30 50 100 200 

Td (ms) 

FIG. 7: Iron mj — J lifetimes in a 1.0 Tesla field vs. diffu- 
sion lifetime at zero field. There is no region of Td for which 
Tj increases. This indicates that Zeeman relaxation is the 
dominant loss mechanism for the entire range of data. 

rules and observe the resulting fit for 7. One extreme 
assumption is Atoj — ±1. The opposite extreme is that 
the rates into all energetically allowed states are equal. 
By fitting the entire data set to the simulation using the 
extremes in selection rules, we obtain the result in Table 
Hin The values of 7 vary between 2 x 10^ and 7 x 10^. 
We assign a lower bound of 7 > 2 x 10'^ based on this 
result. 



B. Iron: Upper Limit on 7 

Fig. [7] shows the measured tj vs. Td- The predicted re- 
gion of Tj cx Td at low buffer-gas density is not observed, 
indicating that Zeeman relaxation occurs on a time scale 
faster than diffusion for the entire range of r^. Therefore 
our finite temperature model cannot be used to fit for 7. 
The slight rise in tj with increasing buffer gas is consis- 
tent with measuring a combination of Zeeman relaxation 
and diffusion as discussed in the previous section. 

We set an upper bound of 7 < 3 x lO'^ by using conser- 
vative values for the parameters in the following expres- 
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7 < TdTRV^G. 



(11) 



We calculate v from the temperature measured via the 
Voigt profile of a zero-field spectrum taken 8 ms after 
ablation. Recall that tj is a reciprocal sum of diffusion 
and Zeeman relaxation lifetimes. When — tji the ex- 
pected value of Tj — Td/2. We assume this condition is 
met at our lowest buffer- gas point since Td/2 is within the 
error bar of tj. This yields the most conservative upper 
bound on 7 because th is at a maximum value consistent 
with our observation that tj < Td- By using the above 
method, we are able to set an upper bound for iron that 
is approximately a factor of 5 smaller than the bound 
measured for scandium ^j. 



C. Impact on Previous Work 

We have demonstrated that when the Zeeman splitting 
between neighboring mj states is comparable to ther- 
mal energy, thermal excitations impact the value of 7 ex- 
tracted from measurements of mj — J decay. First, Zee- 
man relaxation occurs faster than the measured mj = J 
decay, so the extracted 7 is lower than it would be as- 
suming no thermal excitations. Second, uncertainty in 
selection rules for relaxation collisions leads to uncer- 
tainty in how much 7 must be lowered. If we define 
f — ^BgjB/kBT, these effects are significant when ^ < 1. 
The nickel measurement was taken at ^ ~ 0.87. Previous 
work with transition metals pLj] and rare earth atoms 
were performed under similar experimental conditions. 
We now discuss the impact of the present work on those 
measurements. 

The reported value of 7 for titanium, a J — 2 transition 
metal, was found by measuring the decay of the mj — 2 
state. It was assumed that the measured exponential 
decay rate of mj — J atoms was equal to the Zeeman 
relaxation rate so long as the Zeeman degrees of freedom 
were not in thermal equilibrium with the translational 
temperature Ttrans- Specifically, the Zeeman tempera- 
ture Tz was defined at a particular field B by equating 
the ratio of the local populations of two levels, mj and 
m'j, to the Boltzmann factor at temperature Tz'- 



N„ 



N„ 



exp 



/ gjUBjm'j ~ mj)B 



ksTz 



(12) 



It was claimed that when Tz >> Itians, the observed 
mj — J decay was equal to the Zeeman relaxation rate. 
Although this is a necessary condition to measure Zee- 
man relaxation, it is not sufficient to avoid thermal effects 
outlined in the present work. The titanium measurement 
was performed at 3.8 T with 1.8 K atoms. This results 
in ^ = 0.94; very similar to our measurement conditions 
for nickel. Therefore, for a given set of selection rules, we 
expect a comparable downward shift in 7 which would 



result in better agreement between experiment [l| and 
theory [33^1 . As the number of Zeeman levels decreases, 
the number of decay channels decreases and the uncer- 
tainty in selection rules diminishes. Because titanium is 
a J = 2 atom, it has five Zeeman levels compared to 
nine for nickel. Therefore, the uncertainty introduced by 
selection rules for titanium will be less than those found 
for nickel. 

The rare earth 7 values were found by measuring the 
decay of a trapped atomic sample. Because multiple low- 
field seeking Zeeman states were simultaneously trapped 
in an inhomogeneous magnetic field, it was impossible to 
monitor the decay of the mj = J state via an isolated 
spectroscopic line. As a result, a model which included 
thermal excitations, atom drift due to the trapping field, 
and diffusion was implemented to simulate the dynamics 
of all Zeeman states. Therefore, the reported 7 values 
already take into account thermal effects. However, un- 
certainties in selection rules were not addressed. It was 
assumed that toj = J atoms could decay into any en- 
ergetically allowed state with equal probability (the sec- 
ond row of Table UTT)) . The values of J for the rare earth 
elements studied range from 7/2 (thulium) to 8 (dyspro- 
sium) compared to 4 for nickel. As a result, uncertainties 
in 7 due to unknown selection rules should be comparable 
to or worse than those found for nickel. 



V. SUMMARY AND CONCLUSIONS 

A study of collisions between the most low-field seek- 
ing Zeeman state of nickel and iron with '^He has been 
performed to determine the feasibility of buffer-gas load- 
ing highly magnetic transition metals into a magnetic 
trap. Atoms were introduced via laser ablation into a 
cryogenic cell containing a background gas of '^He. Al- 
though we could not measure buffer-gas density, its rel- 
ative density was measured by observing the diffusion 
rate of the atomic sample through the buffer gas. We 
measured the ratio 7 of diffusion cross-section to angular 
momentum reorientation cross-section by measuring the 
mj = J state decay at several different buffer-gas den- 
sities. For our experimental conditions the energy split- 
ting between adjacent Zeeman levels was comparable to 
the thermal energy of the atomic sample. Our operating 
temperature was set by the ablation power required to 
achieve an adequate signal-to-noise ratio. We could not 
operate at higher fields because magnetic broadening of 
the atomic resonances decreased our signal-to-noise to 
intolerable levels. Under these circumstances, thermal 
excitations into the nij — J state cause its decay to dif- 
fer from pure Zeeman relaxation. In order to find an 
accurate value of 7, we fit measured mj = J decay to a 
model of Zeeman state dynamics that includes thermal 
excitations. 

For the Ni-'^He system 7 was found to be Sl^'g (stat) ± 
l{sys) X 10'^, assuming the Zeeman state relaxation cou- 
pling coefficients of Table HTl A change in the assumed re- 
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laxation coupling coefBcients changes the predicted ther- 
mal excitation rates into the mj = J state, resulting 
in further uncertainty in 7 as listed in Table IIIII We 
have also set an upper limit on 7 for the Fe-'^He system 
of 3 X 10'^. These values of 7 are high enough to allow 
buffer-gas cooling to thermalize Ni and Fe, but too small 
to allow a sizeable sample to remain trapped after the 
100 ms required for removal of the buffer gas. 

Our 7 measurements extend the experimentally ex- 
plored range of transition metal-helium collisions to 
species with high magnetic moments. Our method of 
finding 7 by measuring toj — J decay at many buffer-gas 
densities allows us to measure smaller values of 7 and set 
smaller upper bounds than in previous work. While we 



do find that inelastic angular momentum changing col- 
lisions are suppressed in the Ni-He and Fe-He systems, 
the degree of suppression is low compared to the rare 
earth elements, which have 7 > 10^ [21 • The values for 
7 measured here are similar to the values measured for 
other transition metals, specifically the Sc-He and Ti- 
He systems ^Ij. Our observations, together with these 
previous measurements, are consistent with a hypothesis 
of reduced collisional angular momentum transfer due 
to screening of the valence electrons by closed electron 
shells. 

This work was supported by the Office of Naval Re- 
search, the National Science Foundation, and the NSF 
Harvard/MIT Center for Ultracold Atoms. 
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